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A new contribution from pair-production of light gluinos, of mass 12 to 16 GeV, with two-body decays into bottom quarks and light 
bottom squarks, helps to obtain a bottom-quark production rate in better agreement with hadron collider data. The masses of the 
gluino and bottom squark are further restricted by the ratio of like-sign to opposite-sign leptons at hadron colliders. Constraints on 
this scenario from other data are examined and predictions are made. 
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1 Motivation and Constraints 

The cross section for bottom-quark production at hadron 
collider energies exceeds the central value of predictions 
of next-to-leading order (NLO) perturbative quantum 
chromodynamics (QCD) by about a factor of twoEl. This 
longstanding discrepancy has resistedJuUy satisfactory 
resolution within the standard model □. The NLO con- 
tributions are large, and it is not excluded that further 
higher-order effects in production and/or fragmentation 
may resolve the discrepancy. However, the disagreement 
is surprising because the relatively large mass of the bot- 
tom quark sets a hard scattering scale at which fixed- 
order perturbative QCD computations of other processes 
are generally successful. The data invite the possibility 
of a contribution from "new physics" . 

In a recent paper i, my collaborators and I explore 
an explanation within the context of the minimal super- 
symmetric standard model. We postulate the existence 
of a relatively light gluino g (mass ~ 12 to 16 GeV) that 
decays into a bottom quark b and a light bottom squark b 
(mass ~ 2 to 5.5 GeV). The g and the b are the spin-1/2 
and spin-0 supersymmetric partners of the gluon (g) and 
bottom quark. In our scenario the masses of all other su- 
persymmetric particles are arbitrarily heavy, i.e., of order 
the electroweak scale or greater. While the gluino decays 
instantly with 100% branching fraction into the b and 
b, the b is either long-lived or decays via R-parity vio- 
lation into a pair of hadronic jets. We obtain improved 
agreement with hadron collider rates of bottom-quark 
production. Several predictions are made that can be 
tested readily with forthcoming data. 

Our assumptions are consistent with all constraints 
on thEpinasses and couplings of supersyrnmetric parti- 
cles □'LIeII3^eI. If the light bottom squark (6i) is an appro- 
priate mixture of left-handed and right-handed bottom 
squarks, its tree-level coupling to the neutral gauge boson 
Z can be small, leading to good agreement with the Z- 
peak observablesQ. The couplings Zi^i^ and Z^^^j^^, where 

62 is the heavier bottom squark, survive but present 



no difficulty if 77ig^ > 200 GeV. Bottom squarks make 
a small contribution to the inclusive cross section for 
e~^e~ hadrons, in comparison to the contributions 
from quark production, and bb resonances axe likely to 
be impossible to extract from backgrounds The an- 
gular distribution of hadronic jets produced in e'^e~ an- 
nihilation can be examined in order to bound the con- 
tribution of scalar-quark production. Spin-1/2 quarks 
and spin-0 squarks emerge with different distributions, 
(1 ± cos^S), respectively. We re-fit the angular distri- 
bution measured by the CELLO collaboration Q, and 
find it is consistent with the production of a single pair 
of charge- 1/3 squarks along with five flavors of quark- 
antiquark pairs. A new examination of the angular dis- 
tribution with greater statistics woulibe valuable. The 
exclusion by the CLEG collaboration^ of a 6 with mass 
3.5 to 4.5 GeV does not apply since that analysis focuses 
on the decays b —^ cW and b — > cl. On the other hand, 
these data might be reinterpreted in terms of a bound on 
the R-parity violating lepton-numbcr violating decay of b 
into cl. It would be interesting to study the hadronic de- 
cays b — !■ eg, with q — d or s, and b ^ us with the GLEO 
data. In summary, measurements at e^e" colliders do 
not significantly constrain b masses. An analysis of 2- and 
4-jet events by the ALEPH collaboration eI disfavors g's 
with mass rrig < 6.3 GeV but not g's in the mass raxige of 
interest to us. A renormalization group argument □ sug- 
gests that the existence of a light b goes hand-in-hand 
with a comparatively light g {rUg ~ 10 GeV). 



2 Differential Cross Section 

Because the excess production rate is observed in all 
bottom-quark decay channels, an explanation in terms of 
new physics is guided towards hypothesized new particles 
that decay either like bottom quarks or directly to bot- 
tom Quarks. The former is difficult to implement success- 
fully □. In our scenario, light gluinos are produced in pairs 
via standard QGD subprocesses, dominantly g+g 9+9 
at Tevatron and Large Hadron Collider (LHG) energies. 
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The g has a strong color couphng to &'s and 6's and, as 



long as its mass satisfies rrig > rrih 



the g decays 



promptly to b + b. The magnitude of the b cross section, 
the shape of the 6's transverse momentum pxb distribu- 
tion, and the CDF measurement t3 of i?'' — mixing 
are three features of the data that help to establish the 
preferred masses of the g and b. 

We include contributions from both q + q 9 + 9 
and 9 + g — > 9 + g- The subprocess g + b 9 + b 
contributes insignificantly. In Fig. 1 we show the in- 
tegrated pTb distribution of the b quarks that results 



from g ^ b + b, 



for rrig 



=14 GeV and mr = 3.5 GeV. 



The results are compared with the cross section obtained 
from next-to-leading order (NLO) perturbative QGD and 
CTEQ4M parton distribution functions (PDF's)ll2l, with 
nib = 4.75 GeV, and a renormalization and factoriza- 
tion scale — \/ml + p^^,- SUSY-QCD corrections to 
bb production are not included as they are not available 
and are generally expected to be somewhat smaller than 
the standard QCD corrections. A fully differential NLO 
calculation of ^-pair production and decay does not exist 
either. Therefore, we compute the ^-pair cross section 
from the leading order (LO) matrix element with NLO 

PDF's0, ^ = y^m? -I- , a two-loop a^, and multiply 

by 1.9, the ratio of inclusive NLO to LO cross sectionsEl. 

A relatively light gluino is necessary in order to ob- 
tain a bottom-quark cross section comparable in mag- 
nitude to the pure QCD component. Values of nig ~ 
12 to 16 GeV are chosen because the resulting g decays 
produce pxb spectra that are enhanced primarily in the 
neighborhood of p™^" — nig where the data show the 
most prominent enhancement above the QCD expecta- 
tion. Larger values of rrig yield too little cross section 
to be of interest, and smaller values produce more cross 
section than seems tolerated by the ratio of like-sign to 
opposite-sign leptons from b decay, as discussed below. 
The choice of has an impact on the kinematics of the 
b. After selections on p^b"; large values of reduce the 
cross section and, in addition, lead to shapes of the pxb 
distribution that agree less well with the data. The val- 



are correlated; similar results to those 

12 GeV, but 



ues of TO^ and 

shown in Fig. 1 can be obtained with rrig 
mi ~ rrib. 

After the contributions of the NLO QCD and SUSY 
components are added, the magnitude of the bottom- 
quark cross section and the shape of the integrated p™^" 
distribution are described well. A theoretical uncertainty 
of roughly ±30% may be assigned to the final solid curve, 
associated with variation of the b mass, the renormal- 
ization/factorization scale, and the parton distributions. 
The SUSY process produces bottom quarks in a four- 
body final state and thus their momentum correlations 
are different from those of QCD. Angular correlations 
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Figure 1: Bottom-quark cross section in pp collisions at ^/S = 1.8 
TeV for pxb > Pxb^ with a gluino of mass nig = 14 GeV and a bot- 
tom squark of mass m^^ = 3.5 GeV. The dashed curve is the central 
value of the NLO QCD prediction. The dotted curve shows the Pt 
spectrum of the b from the supersymmetry (SUSY) processes. The 
solid curve is the sum of the QCD and SUSY components. Data 
are from Ref. 1. 



between rcuipiis that arise from decays of 6's have been 
measured E^tj. Examining the angular correlations be- 
tween 6's in the SUSY case, we find they are nearly indis- 
tinguishable from those of QCD once experimental cuts 
are applied. 



3 Same-sign to Opposite-sign Leptons 

Since the g is a Majorana particle, its decay yields both 
quarks and antiquarks. Gluino pair production and sub- 
sequent decay to 6's will generate 66 and 66 pairs, as 
well as the 66 final states that appear in QCD produc- 
tion. When a gluino is highly relativistic, its helicity 
is nearly the same as its chirality. Therefore, selection 
of ^'s whose transverse momentum is greater than their 
mass will reduce the number of like-sign 6's. In the in- 
termediate pt region, however, the like-sign suppression 
is reduced. The cuts chosen in current hadron collider 
experiments for measurement of the ratio of like-sign to 
opposite-sign muons result in primarily unpolarized g's, 
and, independent of the 6 mixing angle, an equal number 
of like-sign and opposite-sign 6's is expected at produc- 
tion. Our SUSY mechanism leads therefore to an increase 
of like-sign leptons in the final state after semi-leptonic 
decays of the 6 and 6 quarks. This increase could be 
confused with an enhanced rate of — & mixing. 
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Time-integrated mixing analyses of lepton pairs ob- 
served at hadron colliders are interpreted in terms of the 
quantity x = fdXd + fsXs, where fd and fs are the frac- 
tions of -B° and hadrons in the sample of semi-leptonic 
B decays, and Xf is the time-integrated mixing probabil- 
ity for B^. Conventional bb pair production determines 
the quantity LSc — 2^(1 — x), the fraction of bb pairs 
that decay into like-sign leptons. Our SUSY mechanism 
leads to a new expression 



LSc 



fqcd 



2 (Tgg -I- fJqcd CTgg-t-CTqcd 



2xcff(l - Xcff). 



(1) 

The factor 1/2 arises because N{bb + bb) ~ N{bb) in the 
new mechanism. Introducing G = (Jgg /cTqcd, the ratio of 
SUSY and QCD bottom-quark cross sections after cuts, 
and solving for the effective mixing parameter, we obtain: 



Xoff 



X 
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VT+G 2 
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(2) 



We interpret the CDF collaboration's measupment as a 
determination of Xoff = 0.131 ± 0.02 ± O.OIGEI and note 
that it is marginall^f larger than the world average value 
X = 0.118 iO.OOSEJ. 

To estimate jmr theoretical Xcfi, we assume that 
X = 0.118±0.005EJ represents the contribution from only 
the QCD bb component. We determine the ratio G in the 
region of phase space where the measurement is niadeE£l, 
with both final &'s having pTb > 6.5 GeV and rapidity 
\yb\ < 1- For gluino masses of nig = 14 and 16 GeV, we 
obtain G = 0.37 and 0.28, respectively, with rn^; = 3.5 
GeV. We compute Xoff = 0.17 for irig = 14 GeV, and 
Xcff = 0.16 with rrig = 16 GeV. To estimate the uncer- 
tainty on G, we vary the renormalization/factorization 
scale at which the cross sections are evaluated between 
fi = Tfixjl and /i = 2TOa;. Uncertainties of ±50% are ob- 
tained and lead to uncertainties in Xcff of i5xoff — ±0.02. 
Additional theoretical uncertainties arise because there is 
no fully differential NLO calculation of gluino production 
and subsequent decay to 6's. 

Comparing our expectations with the CDF value, we 
conclude that values of m-g > 12 GeV lead to a calculated 
Xoff that is consistent with the data within experimental 
and theoretical uncertainties. With Ogg/a^cd ~ 1/3, we 
can satisfy the magnitude and px dependence of the b 
production cross section and the mixing data. 



4 Implications and Remarks 

^.1 Hadron Reactions 

Among the predictions of this SUSY scenario, the most 
clearcut is pair production of like-sign charged B mesons 
at hadron colliders, B^B^ and B~B~. To verify the 



underlying premise of this work, that the cross section 
exceeds expectations of conventional perturbative QCD, 
a new measurement of the absolute rate for b production 
in run II of the Tevatron is important. A very precise 
measurement of x in run II is obviously desirable. Since 
the fraction of b's from gluinos changes with pTb, we also 
expect a change of x with the cut on pxb- The b jet from 
g decay into bb will contain the b, implying unusual mate- 
rial associated with the b in some fraction of the bb data 
sample. The existence of light 6's means that they will be 
pair-produced in partonic processes, leading to a slight 
increase (~ 1%) in the hadronic dijet rate. Our approach 
increases the b production rate at HERA and in 77 col- 
lisions at LEP by a small amount, not enough perhaps 
if earljj-epqperimental indications in these cases are con- 
firmed Ejlla, but a full NLO study should be undertaken 
and better parton densities of photons are needed. 

4-2 Running of as 

In the standard model, a global fit to all observables 
provides an indirect measurement of the strong coupling 
strength as at the scale of the Z boson mass Mz- The 
value asLMz) — 0.119±0.006 describes most observables 
properlyEZl. A light g with mass about 15 GeV and a light 
b slow the running of and modify as{Mz), determined 
by extrapolation from experiments performed at energies 
lower than rUg. The result is a shift of 0.007 in as derived 
from these experiments, to as{Mz) — 0.125, at the up- 
per edge of the uncertainty band. Greater precision on 
the determination of as{Mz) from evolution would be 
valuable. Slower running of as{Q) also means a slower 
evolution of parton densities at small x, an effect that 
might be seen in HERA data for Q > nig. Presence of 
a scalar b in the proton breaks the Callan-Gross relation 
and yields a non-zero leading-twist longitudinal structure 
function Fl(x,Q) at leading- order. 



4.3 b-onia 

Possible bound states of bottom squark pairs could be 
seen as ~ 0^, 1~, 2^, ... mesonic resonances in 77 
reactions and in pp formation, with masses in the 4 to 
10 GeV range. They could show up as narrow states 
in the fJ-^ fJ.~ invariant mass spectra at hadron colliders, 
between the J/'^ and T. At an e+e~ collider, the in- 
termediate photon requires production of a J^^ = 1 
state. Bound states of low mass squarks with charge 2/3 
were studied with a potential modelEl. The small leptonic 
widths were found to preclude bounds for m~ > 3 GeV. 
For bottom squarks with charge —1/3, the situation is 
more difficult. 
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4-4 b lifetime and observability 

Strict R-parity conservation in the MSSM forbids b de- 
cay unless there is a hghter supersymmetric particle. 
R-parity-violating and lepton-number-violating decay of 
the 6into at least one lepton is disfavored by the CLEO 
data Q and would imply the presence of an extra lep- 
ton, albeit soft, in some fraction of b jets observed 
at hadron colliders. The baryon-number-violating R- 
parity- violating (ftp) term in the MSSM superpotential 
is = KikUfD^pi; Uf and £)f are right-handed- 

quark singlet chiral superfields; and fc are genera- 
tion indices. The limits on individual ftp and baryon- 
number violating couplinga_.A" are relatively weak for 
third-generation squarksEaO, A-'^j, < 0.5 to 1. 

The possible ftp decay channels for the b are 123 : 
6 ^ u + s; 213 : c + d; and 223 : 5 ^ c + s. The 
hadronic width is 



T{b 



• jet jet) = -^sin 



(3) 



If mi = 3.5 GeV, r(6 ^ ij) 
all A'' 



0.08|A^'^3|2 GeV. Unless 

•Vy 3 are extremely small, the b will decay quickly and 
leave soft jets in the cone around the b. 5-jets with an 
extra c are possibly disfavored by CDF, but a detailed 
simulation is needed. 

If the b is relatively stable, the b could pick up a light 
u or d and become a or "mesino" with J = 1/2, 
the superpartncr of the B meson. The mass of the mesino 
would fall roughly in the range 3 to 7 GeV for the interval 
of b masses we consider. The charged mesino could fake 
a heavy muon if its hadronic cross section is small and if 
it survives passage through the hadron calorimeter and 
exits the muon chambers. Extra muon-like tracks would 
then appear in a fraction of the bb event sample, but 
tracks that left some activity in the hadron calorimeter. 
The mesino has baryon number zero but acts like a heavy 
p - perhaps detectable with a time-of-flight apparatus. A 
long-lived b is not excluded by conventionaLsearches at 
hadron and lepton colliders, but an analysis o similar to 
that for g's should be done to verify that there are no 
additional constraints on the allowed range of b masses 
and lifetimes. 
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